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Abstract 

Three decades ago Heath found the integral form of the exact an- 
alytic growing mode solution of the linear density perturbation 5 on 
sub-horizon scales including the cosmological constant or the curvature 
term. Recently, we obtained the exact analytic form of this solution in 
our previous work [1] . Interestingly, we are able to extend this solution 
for general dark energy models with the constant equation of state LOde 
in a flat universe. This analytic solution provides the accurate and effi- 
cient tools for probing the properties of dark energy models such as the 
behavior of the growth factor and the growth index. We investigate the 
growth index and its parameter at any epoch with this exact solution 
for different dark energy models and find that the growth index is quite 
model dependent in the redshift space, 0.25 < z < 1.5, so observations 
of the structure growth around this epoch would be very interesting. 
Also one may be able to rule out some dark energy models by using 
the analysis from this exact solution. Thus, the analytic solution for 
the growth factor provides the very useful tools for future observations 
to constrain the exact values of observational quantities at any epoch 
related to the growth factor in the dark energy models. 
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The background evolution equations in a flat Friedmann-Robertson- 



Walker universe {pm + Pde = Per) are 
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where ujde is the equation of state (eos) of dark energy, pm and p^e are the 
energy densities of the matter and the dark energy, respectively. We consider 
the constant w^e- The sub-horizon scales linear perturbation equation with 
respect to the scale factor a is well known [2], given by 
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We are able to find the exact analytic growing mode solution of 5 for any 
value of the constant uj^e- After replacing new parameters Y = Qa^^'^^ and 
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in Eq. (3), we get 
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A trial solution is 5{Y) = Y'^BiY) because it is the most general combina- 
tion of the solution for the above equation (4) . We replace 5 into Eq (4) to 

get 
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The above equation becomes the so called "hyper geometric" equation when 
we replace Y = —Y with the complete solution [3], 
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where F is the hypergeometric function. 
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Figure 1: Behavior of ^ as a function of ujde when = 0.3 and = 0.1. 



Thus, the exact analytic solution of the sub-horizon scales linear perturba- 
tion becomes 

6(Y) 11^13 1 

^ 3a;de ' 2wde ' 2 6a;rfe ' ^ ^ 

This analytic solution does not have any physical meaning before we fix 
the coefficients ci and C2- If we want to have the correct growing mode 
solution from the above analytic solution, then this solution should follow 
the behavior of growing mode solution at an early epoch, say Oj ~ 0.1. In 
other words, the coefficients of the general solution should be fixed by using 
the initial conditions of the growth factor. 
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For example, (ci,C2) = (-0.716894,1.07822) for uJde = -1-2. We do need 
these accurate numbers to show the proper growing mode behavior. After 
we fix the coefficients from the initial conditions, we are able to determine 
the growth factor 5g from the general form of solution 5. 
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In the literature, it is used to normalizing the growth mode with respect 
to its present value. We can obtain from Eq. (7), where lo(= Q) 

means the value of Y at the present time. Therefore, we only need to know 
the ratio of ^ to determine the growth mode with the initial conditions Eq. 
(8). From Eqs. (7) and (8), it is straightforward to show that 
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Even though 3^ is a function of three quantities Oj, Q {i.e. J^^), and tOde-, its 
dependence on ai and Q is quite small. When we vary both ai and from 
0.01 to 0.1 and from 0.2 to 0.35, respectively, the changes on ^ are about 
10~^ % both cases. In Fig. 1, the behavior of the ratio of two coefficients in 
the growth mode solution is depicted as a function of Ude- As ojde decreases, 
the magnitude of the ratio ^ also decreases. This can be better understood 

^ C2 

with Fig. 2. 



In Fig. 2, we show the behaviors of the growth factor 5g for the different 
dark energy models (i.e. for the different values of oJde) when = 0.3. 
The dotted, solid, and dashed lines correspond to Ude = —1.2,-1.0 and 
—0.9, respectively. As Ude decreases, 5g maintains the linear growth factor 
proportional to a for a longer time. This is easy to understand. If there is 
the same amount of the matter at the present epoch for the different models, 
then there will be more matter component in the past for the smaller value 
of OJde- Thus, dark energy model with the smaller value of ujde maintains the 
longer linear growth behavior. 

We are also able to obtain the exact value of the growth index /(a) = 
in any epoch. We also get the exact value of the growth index param- 
eter 7(a) = ijijf^^- We investigate the behaviors of /(a) and 7(a) based 
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Figure 2: Evolutions of 5g{a) when = 0.3 for ojde = —1-2, —1.0, and —0.9 
(from top to bottom). 

on this exact analytic growth factor for the different dark energy models as 
shown in Fig. 2. 

In Fig. 3a, we show the cosmological evolution of f{a) for the differ- 
ent dark energy models when = 0.3. Dashed, solid, and dotted lines 
correspond to oJde = —0.9, —1.0 and —1.2 dark energy models, respectively. 
As we show in Fig. 3a, the smaller value of uj^e gives the larger value of 
the present growth index /(a = 1). The present values of the growth in- 
dex are 0.511,0.513, and 0.516 for corresponding models, respectively. Thus 
we may not see any differences between different models from the present 
growth index. If we investigate the / values at z ~ 0.15 (i.e. a ~ 0.87), 
then the growth index values will be f{z = 0.15) = 0.588,0.598, and 0.618. 
The current 2dF observation value is 0.36 < f{z = 0.15) < 0.66 [4, 5]. 



In Fig. 3a, the behavior of /(a) shows stronger model dependence around 
0.25 < z < 1.5. Thus, future observations in this redshift range should 
give us better observational constraints on the growth index. In Fig. 3b, 
we show the evolution of the growth index parameter 7(a) for correspond- 
ing models with the same notation as Fig. 3a. The present values 7(0 = 
1) = 0.5581,0.5547, and 0.5491 for ujde = -0.9,-1.0, and -1.2, respec- 
tively. There is only 2% deviation from the different models. We are 
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Figure 3: a) Evolution of /(a) when fij^ = 0.3. b) 7(a) evolution with the 
same VL^. Dashed, solid, and dotted lines correspond to ujde = —0.9,-1.0, 
and —1.2, respectively. 

able to obtain the 7 values at 2: = 0.15 (equally a ~ 0.87). We get 
7(2: = 0.15) = 0.5559,0.5525, and 0.5468 for corresponding models. This 
shows not much changing in 7 between the models. However, 7 values at 
high z show very strong dependence and vary dramatically [1]. Thus, 
we had better to limit 2; < 2 ~ 3 for the investigation of 7. We summarize 
the results in Table 1. We also indicate the coefficients of the growth factor 
in this table, which are obtained from the growth factor initial conditions 
Eq. (8). In the last column, we compare the present growth index param- 
eter 7{y5 obtained from Ref. [6]. Those values are quite close to the ones 
obtained from the exact analytic solution. Thus, the approximation of 7 
given in Ref. [6] is a good one as long as one probes 7 near the present. 



In addition to this, if we naively take the surface values of the l-cr result 
of the 2dF measurement without taking into account the uncertainties in 
the relevant cosmological parameters, then we may be able to rule out many 
dark energy models by using this exact analytic solution. For example, if 
LOde = —3 like in some phantom models, we will obtain f{z = 0.15) = 0.766 
for = 0.3, which is way too large compared to the 2dF observational 
value even when we consider the error in the data coming from the selection 
effects. Thus, dark energy models with too small values of uJde should be 
ruled out if we want to keep the concordance model. However, there are 
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Table 1: ojde is the eos of the dark energy. 7(^=0) f{z=o.i5) ^re the 
values of the growth index at the present and z = 0.15, respectively. 7(^=0) 
and 7(2=0.15) sue growth index parameters at the corresponding epoches. ci 
and C2 are the coefficients of the growth factor 6g obtained from the initial 
conditions in Eq. (8). 7^5 are the values of the growth index parameters 
obtained from Ref. [6] . 

still large uncertainties in the 2dF measurement and we need to wait for 
more accurate measurements to conclude this. Also some models with time 
varying code might have better chance to survive. 

It has been misunderstood that the sub-horizon scale growth factor for 
the general constant ujde is known [7]. However, the solution in the given 
article is not correct because it claims that the growth and decaying solution 
are separable for the general oJde- It is well known that the separation of 
two modes is possible only when ujde = ~5 or — 1 [8]. The detail discussion 
about this is out of the main stream of letter [9] . 

Even though the results of this letter is limited for the constant equation 
of state of dark energy, we are able to apply these solutions to the time- 
varying oJde by interpolating between models with constant uJde [9]. Also 
it is well known that the time-dependence of oJde is extremely difficult to 
discern because the dark energy is dynamically unimportant at the redshifts 
where oJde departs from its low z value. In addition, for the substantial 
changes in oJde at low rcdshift, there is always a constant LOde that produces 
very similar evolution of all of the observables simultaneously [10, 11]. Also 
this analytic solution can provide useful templates to study the structure 
growth in dark energy models with time varying equation of state. 

In our previous work, we also obtained the exact analytic solution of 
the growth factor 5g for the cosmological constant case [1]. Even though, 
dg seems to be quite different from the solution dg in this letter, both are 
indeed same solutions which show the same physical behaviors. The details 
of comparison of them are irrelevant to the results of this letter and we do 
not show any detail in the present consideration [9, 12]. 
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